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Abstract 

A multidimensional chaos is generated by a special initial value prob- 
lem for the non-autonomous impulsive differential equation. The existence 
of a chaotic attractor is shown, where density of periodic solutions, sensi- 
tivity of solutions and existence of a trajectory dense in the set of all orbits 
are observed. The chaotic properties of all solutions are discussed. An ap- 
propriate example is constructed, where the intermittency phenomenon is 
indicated. The results of the paper are illustrating that impulsive differen- 
tial equations may play a special role in the investigation of the complex 
behavior of dynamical systems, different from that played by continuous 
dynamics. 

1 Introduction and Preliminaries 

The investigation of the complex behavior of dynamical systems, as well as the 
development of the methods of this investigation, has made a long way, starting 
with significant works [1]- [T5] . 

It is natural to discover a chaos [l]-[3],[6], [7]- [12], [El [IS], and proceed by 
producing basic definitions and creating the theory. On the other hand, one 
can, following the prescriptions, shape an irregular process by inserting chaotic 
elements in a system which has regular dynamics otherwise (let's say is asymp- 
totically stable, has a global attractor, etc). This approach to the problem also 
deserves consideration as it may allow for a more rigorous treatment of the phe- 
nomenon, and help develop new methods of investigation. Our results are of 
this type. Using the logistic map as a generator of moments of impulses in the 
multidimensional system, we observe the density of the periodic solutions, the 
sensitivity of solutions and the existence of a trajectory dense in the set of all 
solutions in a bounded region of the space. 

The first mathematical definition of chaos was introduced by Li and Yorke 
[9]. They proved that if a map on an interval had a point of period three, then it 
had points of all periods. Moreover, there exists an uncountable scrambled subset 
of the interval. Devaney [^ gave an explicit definition of a chaotic invariant set 
in an attempt to clarify the notion of chaos. To the properties of transitivity and 
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sensitivity [12] he added the assumption that the periodic points are dense in the 
space. An intensive discussion of the definitions has continued during the last 
decades. It was shown in [16] that a map is Li-Yorke chaotic if and only if there 
is a two-point scrambled set. Paper [17] proves that chaos as defined by Devaney 
is stronger than that defined by Li-Yorke. In [18] it was shown that transitivity 
and density of periodic points imply sensitivity. 

In this paper we concentrate on the topological components of the version 
proposed by Devaney. The special initial value problem is introduced, when 
the moments of the impulsive action are functionally dependent on the initial 
moment. 

One of the most powerful tools of the chaos investigation is the conjugacy 
with the symbolic dynamics. Our results are based on the method, too. Consider 
the sequence space [4] 



where t = (tg^i . . .) G S2, and the shift map a : II2 —>■ S2, such that a{s) = 
{siS2 . . .). The semidynamics (S2, a) is the symbolic dynamics |19j . 

The map is continuous, cardPern{cr) = 2"-,Per{a) is dense in E2, and there 
exists a dense orbit in S2. 

We assume that there exist a homeomorphism S between S2 and a set A C 
1,1 = [0,u;], where is a fixed positive number, and a map h : A ^ A, such 
that S o h = a o S. That is h and a are topologically conjugate. It is known 
that S2 and A are Cantor sets: they are closed, perfect and totally disconnected 
[19] . Obviously, they are compact. Moreover, a and h have sensitive dependence 
on the initial conditions, periodic points are dense in S2 and A respectively, and 



S2 = {s = (S0S1S2 . . .) : Sj = or 1} 



with the metric 
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the maps are topologically transitive. We may also specify that for each p G N 
there exists a solution with period p, and topological transitivity means having 
a dense trajectory. 

One of the most popular examples of the map h is the logistic map /ia;^(l — 
x),yU > 4, considered on a subset of [0, 1], [20]. 

The description of the main subject of our paper should begin with the dis- 
cussion of the moments of impulses, as their generation is most important for 
the emergence of chaos. 

For every E A one can construct a sequence /t(to) of real numbers Kj, z G Z, 
in the following way. If z > 0, then Kj+i = h{Ki) and = to- Let us show, how 
the sequence is defined for negative i. Denote s° = Sito), s° = (sqS^ . . .). Consider 
elements s = {Os^sl . . .),s = (IsgsJ . . .) of E2, such that a{s) = a{s) = s° and 
t = S^^{s),t = S'^^(s). The homeomorphism implies that hit) = h{t) = to- Set 
h'^ito) may consist of not more than two elements t, t G A. Each of these two 
values can be chosen as K_i(to)- Obviously, one can continue the process to —00, 
choosing always one element from the set h~^. We have finalized the construction 
of the sequence, and, moreover, it is proved that /t(to) C A. Thus, infinitely many 
sequences fi:(to) can be constructed for a given t^. However, each of this type of 
sequence is unique for an increasing i. Fix one of the sequences and define a 
sequence C(to) = = + G Z. The sequence has a periodicity property 

if there exists p G N such that Qj^p = + P^i ^ ^- If 'we denote by 11 the set 
of all such sequences {Ci\-,i ^ Z, then a multivalued functional w : I ^ Tl such 
that Ciio) = ""^(^o) is defined. 

In what follows, we assume, without loss of generality, that uj = 1. 

Let J C M be an open interval. We introduce the distance ||C(to) ^C(ti)l|j = 

sup ICi(to) — Ci(ti)l) and, if C(t) = {Ciif)} is a sequence from 11, and m is 

0(to),0(ti)e-/ 

an integer, we denote ((t,m) = {(i+m(t)},i ^ P^- More information about the 
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opportunities connected with this distance can be found in |21j . 

Next, let us consider some useful properties of the elements of 11. They are 
simple consequences of the topological conjugacy. The first one is a reformulation 
of the known property for symbolic dynamics ^1^ . 

Lemma 1.1 Assume that h is topologically conjugate to a then the following 
assertions are valid: 

(a) for each ((to) G 11, an arbitrary small e > 0, and an arbitrary large positive 
number E there exists a sequence C{ti) G 11 with the periodicity property 
such that IC(io) — C('^i)|j < where J = {0,E) 

(b) There exists a sequence ({t*) G 11 such that for each to G A, and an arbi- 
trary small e > 0, and an arbitrary large positive number E there exists an 
integer m such that \((to) — ((t*,m)\j < e, where J = {0,E). 

Proof: (a) Fix numbers tg & A,e > 0, E > 0, and denote S'(to) = s'^ = 
{sqSi . . .). Since S" is a homeomorphism and E2 is compact, there exists a number 
6>0 such that \S-\s^) - S-\s^)\ < e if d[s\ s^] < 5, where s\ G Eg. 

Next we define a periodic sequence s = (sqSiSs . . .) from E2 to satisfy e and 
E. Take a number / G N such that / > E and < Assume that Sj = i = 
0, 1, . . . , 2/ — 1, and = Si,i & Z. Consider the sequence Kj(to) = h'^{to),i > 0. 
Since h\to) = S^^ o a* o S", and d[a^s^,a^s] < 6,i = 0,1,...,/, we have that 
|/i*(to) - h''{S''^{s)\ < e. So, if we denote t = S~'^{s) then |Kj(to) - K.i{t)\ < e,i = 

0. 1. ... ,1. In other words, ||C(^o) — C(^)ll J < where J = (0, E). The assertion 
is proved. 

(6). Consider the sequence s* G E2 such that 

s* = ^ I poonoii I . . . , 

1 element blocks 2 element blocks 

that is s* is constructed by successively listing all blocks of O's and Vs of length 

1, then length 2, etc. Sequence s* is dense in E2 [1]. Denote t* = S~^{s*). Let 
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us fix to G A, e > and E > 0. Similarly to the previous proof, fix a number 
5>0 such that \S-\s^) - S-\s^)\ < e if d[s\s'^] < 5, where s\s'^ e E2. Take 
I e N such that I > E + 1 and < 5. We can find 21— block s^s^.^^ . . . s'^_^_2i 
such that i— th element of S(to),i = 0, 1, ... ,2/, is equal to s*_^^. One can see 
that IIC(^o) - C(^*)"^)IU < ^> where J = (0, E). The lemma is proved. 

The following special initial value problem for the impulsive differential equa- 
tion, 

z'{t)=Az{t)+f{z),tj^ato), 
A\t=c,^to)-Bz{Ci{to)) + W{z{Ci{to))), 

z{to)= Zo, {to:Zo) eAxW, (1) 

where z G M",t G R, and C(^o) = w{to), is the object, which will be mainly 
discussed in our paper. 

We shall need the following basic assumptions for the system: 

(CI) A, B are n x n constant real valued matrices, det(X-|- B) ^ 0, where Z is 
n X n identical matrix; 

(C2) for all xi,X2e M" functions f{x, y) -.W xW ^ R", W -.W ^ R", satisfy 

- /(a;2)|| + \\W{x,) - W{x2)\\ < L\\xi - X2\\, (2) 
where L > is a constant; 

(C3) 

Bx + W{x) y^Oyx eW; (3) 

(C4) functions / and W are uniformly bounded so that 

sup ||/(a;)|| + sup \\W{x)\\ = Mq < 00. (4) 
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We considered the systems of form ([T]) in ^22j in the investigation of blood 
pressure distribution. 

In the present paper the sensitivity, transitivity and existence of dense pe- 
riodic solutions of nonlinear problem ([T]) are considered. We also prove the 
existence of the chaotic attractor, which consists of all solutions bounded on the 
real axis and discuss the period-doubling cascade for the problem. 

All the results are rigorously verified using the results of the theory of impul- 
sive differential equations [23]- [31]. 

While many theorems of the theory of impulsive differential equations are 
obtained by analogy with the assertions which already exist for the systems 
with continuous solutions, the results of the present work do not have direct 
predecessors from the continuous dynamics. 

The investigation is inspired by the discontinuous dynamics of the neural 
information processing in the brain, information communication, and population 
dynamics [32]-[38]. While there are many interesting papers concerned with the 
complex behavior generated by impulses, the rigorous theory of chaotic impulsive 
systems remains far from being complete. Our goal is to develop further the 
theoretical foundations of this area of research. 

The paper is organized in the following manner. The existence of the attractor 
is under discussion in the next section. In section 3 we consider the main subject 
of the paper: ingredients of chaos, the chaotic attractor, the period-doubling 
cascade and an example. 

2 An attractor consisting of the solutions 
bounded on M 

In this section we discuss the existence of bounded and periodic solutions of the 
initial value problem, and the existence of an attractor consisting of the bounded 
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on M solutions. 

Let us start with a system of a more general form than ([1]). 
Fix a sequence 6 = {9i}, i E Z, such that i < 9i < i + l,i E Z, and consider 
the following initial value problem 

z'{t) = Az{t) + f{z),tj^e„ 
A\t=e, = Bz{e,) + Wiz{ei)), 

z{eo) = zo,zoeR\ (5) 

where matrices A and B, functions / and W are the same as in ([T]). 

In order to define the solutions of the impulsive systems we need the following 
spaces of discontinuous functions. 

A left continuous function z{t) is from a set of functions VC^{9, M), if: 

(i) it has discontinuities only at points 9i,i G Z, and they are of the first kind; 

(ii) the derivative z'(t) exists at each point t G ]R\{6'j}, and at points 9i, z G Z, 
the left-sided derivative exists. 

A function z{t) G VC^O,]^) is a solution of if: 

(1) the differential equation is satisfied for z{t) on ]R\{^^j}, and it holds for the 
left derivative of z{t) at every point 6i,i G Z. 

(2) the jumps equation is satisfied by z{t) for every i G Z. 

It is known [25] that under the above mentioned conditions the solution of 
(IS]) exists, is unique on M, and it is an element of VC^{9, M). 
We shall also need the following set of functions. 

A left continuous function z{t) : [Oq, oo) — M" is from a set of functions 
PC\e, [^o,oo)),if: 
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(i) it has discontinuities only at points 9i,i > 0, and these discontinuities are 
of the first kind; 

(ii) the derivative z'{t) exists at each point t G [6Q,oo)\{6i}, and at points 
9i,i > 0, the left-sided derivative exists; 

A solution z{t) of ([5]) on [^O; oo) is a function z(t) E VC^{9, [Oq, oo)) such that: 

(1) the differential equation is satisfied for z(t) on [^^O; oo)\{6'j}, and it holds 
for the left derivative of z{t) at every point 6i > 0. 

(2) the jumps equation is satisfied by z(t) for every i > 0. 

Let us denote Z{t, s) the transition matrix of the associated with linear 
homogeneous system 

z'{t) = Az{t),t^ei 

Az\t=e, = Bz{e,). (6) 

We may assume that: 

(C5) the matrices A and B commute, the real parts of all eigenvalues of the 
matrix A + ln(X + B) are negative. 

Condition (C5) implies, [23], Theorem 34, that there exist two positive numbers 
and uj, which do not depend on 6, 

\\Z{t,s)\\<Ne-'^^'^'\t>s. (7) 

Obviously, that condition for the matrices to commute is strong. One can 
obtain the exponential decay of the transition matrix using other conditions on 
matrices A and B, but this would lengthen the paper significantly, and, moreover, 
decrease the clarity of the proof. 

Assume additionally that 
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(C6) NL[l + j^]<l. 

The proof of following two Theorems 12.11 and 12.21 replicate those of Theorems 37 
and 89 from [25], respectively. 

Theorem 2.1 // conditions (CI), (C2), (C4) - (C6) are valid, then: 

1. there exists a unique hounded solution z{t) of ^ from VC^{9,M), 

and \\z{t)\\ < NMo[^ + jz^] for all t E M; 

2. if there exists a number p E N, such that Oi^p = 9i + p,i E Z, then the 
bounded solution z{t) has the period p; 

We may assume that 

(C7) -U + NL + ln(l + NL) < 0. 

Let z{t) be the bounded solution, which exists by Theorem 12.11 and zi{t) be 
another solution of ([5]). The following theorem is valid. 

Theorem 2.2 Assume that conditions (CI), (C2), (C4) - (C7) are fulfilled. 
Then, for each t, s E M,t > s, 

Mt) - zi{t)\\ < \\z{s) - zi{s)\\{l + ArL)e(-"+^^+^°(^+^^))(*-^). (8) 

That is, the bounded solution z{t) attracts all solutions of 

In what follows we denote z{t,^,v),C, E M,v E M", a solution of ([T]) with 
to = ^,zq = V, and assume that all conditions (CI) — (C7) are fulfilled. 

Assume that in ([5]) 6 = C(^o);^o £ A. Then Theorem 12.11 implies that for 
each to E A there exists a unique bounded on M solution of ([T]). We denote the 
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solution by z{t, to). Moreover, if C(^o) has the periodicity property with period p, 
then the solution z{t, to) is p— periodic. Denote the periodic solutions as to)- 
By Theorem I2.2l the bounded solution z{t, to) attracts all solutions of ([1]) with 
the same points of discontinuity. That is, if z{t,to, zq) is another solution of ([I]), 
then 

||^(t,to)-^(t,to,^o)|| < lk(to,to)-^o||(l + iVL)e(-+^^+i°(i+^^«(*-*°), (9) 

for all t > to. Finally, we remark that there exist infinitely many periodic so- 
lutions 0(t,to),to G A. More precisely, for each p E N there exists to G A such 
that z(t, to) is a p— periodic solution. These periodic solutions are different for 
diferent p since by condition (C3) they have sequences of discontinuity points 
which do not intersect. 

Denote PC = {z{t, to, zq) : to E A, zq e M"}, and let VCB C PC be the subset 
of all solutions bounded on M. That is, VCB = {z{t,tQ) : to G A}. On the basis 
of above made discussion we may say that the bounded set VCB is an attractor 
with the basin VC. 

3 The chaos 

In this section we introduce the topological ingredients of the chaos for the dy- 
namics of ([1]) and find the conditions for the existence of these ingredients. More- 
over, we will prove that VCB is the chaotic attractor, and discuss the period- 
doubling cascade for the problem. 

3.1 The chaotic attractor existence 

Definition 3.1 We say that (C]) is sensitive on VCB if there exist positive real 
numbers eo, ei such that for each to G A, and for every 6 > one could find 
a number ti G A, |to — ti| < 6, and an interval Q from [to, oo) with length no 
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less than ei such that \ \z{t,to) — z(t,ti)\ \ > eo,t G Q, and there are no points of 
discontinuity of z(t,tQ), z(t,ti) in Q. 

To introduce following definitions of density we need the concept of "close- 
ness" for piecewise continuous functions. Different types of metrics and topolo- 
gies for discontinuous functions are described in [23], [25]- [31], [21], [39], [40]. In 
what follows we are going to use the concept from |39j . 

Let us fix to, ti E I and an interval J C [to, oo)n[i(:i, oo), We say that a function 
^(t) G 'PC\Cito),^ato)) is e-equivalent to a function ^/'(t) G ■pC^{Cih),'^ati)) 
on J and write ^{t){e, J)^/j{t) if ||C(to) - C(^i)l|j < e and \\^{t) - ip{t)\\ < e for 
all t from J such that t ^ ^Ci{to),Ci{ti)ej[(i{'to) , where [a, b],a,b E M, means 

an oriented interval, that is [a, b] = [a, 6] if a < b, and [a, b] = [b, a], if otherwise. 
The equivalence of two piecewise continuous functions when e is small means, 
roughly speaking, that they have close discontinuity points, and the values of the 
functions are close at points that do not lie on intervals between the neighbor 
discontinuity points of these functions. 

Definition 3.2 The set of all periodic solutions (j){t) = (j){t, to), to G A, o/ HP is 
called dense in VCB if for every solution z{t,ti),ti G A, and each e > 0, E > 0, 
there exists a periodic solution (f){t,t*),t* G A, and an interval J C with 
length E such that (j){t){e, J)z{t). 

Definition 3.3 A hounded solution z{t) = z{t,t*),t* g A, of ^ is called dense 
in the set of all orbits of VCB if for every solution Zi{t) = z(t,ti),ti G A, of 
(QP, and each e > 0, E > 0, there exists an interval J with length E and a real 
number ^ such that z{t + ^)(e, J)zi(t). 

Definition 3.4 The attractorVCB is chaotic if 

1. (|7]j is sensitive in VCB; 
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2,. the set of all periodic solutions (j){t,to),to G A, is dense in VCB; 
3. there exists a solution z(t,to),to G A, which is dense in VCB. 

Theorem 3.1 Assume that conditions (CI) — (C7) are fulfilled. Then the set 
of all periodic solutions (j){t,tQ),tQ G A, of IjJ^ is dense in VCB. 

Proof: Let us fix ti G A and E,e > 0. By Lemma 11.11 (a) for an arbitrary 
large number T there exists a sequence C(^o) ^ n with periodic property such 
that llC(ti) -C(^o)||q < e, where Q = {ti,ti + f + E). We shall find the number 
T sufficiently large such that the bounded solution z(t) = z(t, ti) is e— equivalent 
to (j){t, to) on J = {ti+f,ti+f + E). Denote Mi = NMo{^ + 

Let also Zi{t,s) = Z{t, s,((ti)) and Z2{t,s) = Z{t, s,Q{to)),t > s, be the 

transition matrices. We have that [25] 

z{t) = Z^it, l)z{l) + / s)f{z{s))ds + E Z^it, 

1 i<0<t 

= z^it, 1)0(1) + / z^it, s)f{<p{s))ds + E ^2(t, ato))w{<p{ato))). 

1 i<Ci<t 
It is difficult to evaluate the difference between z{t) and (f){t) using the last 

two expressions since the moments of discontinuities do not coincide. For this 

reason let us apply the method of S— equivalence developed in papers [26]-[3T]. 

Let us consider the following system 

v'{t) = Av{t) + f{v),t^ato), 
and the system 

z'{t)=Az{t) + f{z),t^Q,(t^), 

^i=at.) = Bz{ati)) + w{at,)), (11) 

where to)^i ai"e the numbers under discussion. Assuming, without loss of gener- 
ality, that Ci(^o) < 0(^i) foi' all J) we introduce 
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Wl{z) = {I + B) (e^K«(*i)-0(*o)) -I)z+ / e^^^^^''^'''^f{z{s))ds + 

where z(t),zi(t) are the solutions of 

z'{t) = Az{t) (12) 

such that z{Q{tQ)) = z and Zi{(i(ti)) = z{(i(ti)+). One can easily verify that 
— ^^P||^||<Mi,iGZ II ^ Systems ([10]) and ([TT]) are 5— equivalent 

[2S]-[2I]- That is their two solutions with the same initial data coincide in their 
common domain if only t ^ {(i{to),Q{ti)],'i G 

So, if v(t),v{l) = z{l), is the solution of (ITUi) . then v{t) = z{t) for all t ^ 

iQ{to)X^ih)],teZ. 

For v{t) we have that 

t 

v{t)=Z2it,l)v{l) + j Z2{t,s)f{v{s))ds+ 

1 

z,{t,at,))[w{v{ah)) + w,{v{ah))]. 

i<Ci<t 

Consequently, 

mt) - vm < 110(1) - v{l)\\\\Z,{t, 1)11 + \\Z,{t, s)\\LU{s) - v{s)\\ds+ 

\\z2{txAm\LU{Cjm-v{cAm\+ E ii^2(t,o(to))iiiii^i(^;(c.(to))ii< 

l<Cj(io)<t l<0(io)<t 

2MiN + M2 ^ [ Ne-'''''"'>L\\z{s) -v{s)\\ds+ 

1 — e ji 

E iVe-(*-^^(*o))L||t;(0(to))-^;(0(to))||. 
i<Cj<t 
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Now, applying the analogue of Gronwall-Bellman Lemma [25] for discontin- 
uous functions, one can find that 



\\z{t) - v{t)\\ < {2MiN + M2 _ _J e^-'-+^L)(t-i) Y[ [1 + NL) < 
{2MiN + M2^^^)e(-"+™'^(i+^^))(*"i) . (13) 

Inequality implies that \\z{t) -v{t)\\ < eift > T,t ^ [Ol^jTSl^i)], ^ e Z, 
where 

f _ I _l_ ^^^2MiAr+M2e"(l-e-"~) 

-uj + NL + \n{l + NL)' 
(we may assume that e < 2Mi ). That is why, z{t){e, J)4>{t) if J = (ti +T,ti + 
T + E). The theorem is proved. 

Theorem 3.2 Assume that conditions (CI) — (C7) are fulfilled. Then there 
exists a solution of (QP z{t,t*) dense in VCB. 

Proof: Fix positive E, e. By Lemma [1.11 (b) there exists t* G A such that C(^*) is 
dense in 11. Denote z*(t) = z{t,t*). Let us prove that 2*(t) is the dense solution. 

Consider an arbitrary solution z{t) = z{t,to) G VCB. Consider an interval 
Ji = (0,i?i), where Ei is an arbitrarily large positive number. By Lemma ll.ll 
ip) there exists a natural m such that 

||C(ti)-C(r,m)||j, <e. (14) 

We have that 

nt+m 

z,{t + m)= e^(*+'"^i-'")^,(l + m) + / e^^'^"'-''^ f{z,{u))du+ 

Jl+m 

e^(*+'"-^'(*«))iy(z,(0(to))) = e^(*-i)^,(l + m) + 

l+m<Ci{to)<t+m 

f e^(*-")/(z.(M + m))du + e^(*+™~^^(*°»l^(z.(0(to))). 

l+m<C,i{ta)<t+m 
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and 



s<Ci(ti)<i 



Now, using the last two formulas, similarly to proof of Theorem 13.11 using 
(fT4ll and the 5— equivalence technique, we can find a sufficiently large number 
El > 2E, and a natural number m such that z^{t + m) and are e— equivalent 
on J = {El/ 2, El). The theorem is proved. 

Denote m = max|u|<i ||e^"||,m = min|„|<i He^^^H. 

Condition (C3) implies that 

rj = min (Bx + W(x)) > 0. 

||i'||<Mi 

Fix a number o > 3, such that - < |=. 

^ — ' q dm 

We shall need the following assumption. 

r 2m 1 -1 

(C8) L < If^. 

^ ' 2Mi (m+m) 

Theorem 3.3 Assume that conditions (CI) — (C8) are fulfilled. Then (QP 
sensitive on VCB. 

Proof: Fix a solution z{t) = z{t, to), to G A, and a positive 5. 

Let S'()f:o) = s° = (50,5?. . .). Fix a number ti G A such that 5'(ti) = S"*^ = 
(sg, . . . , 4, 4+1' ■5°+2' • • •)> 4 7^ 4' for some > 0. We have that 



^ ' J I lii > n. 



Assume that n > 3, n is sufficiently large for 1^0—^1 1 = — S'^-'^(s-'^)| < 6. 

Now, denote Zi{t) = z(t,ti) the solution of ([1]). 

Since 5* is a homeomorphism and set S2 is compact, for a given i,0 < i < n, 
the set 

P, = {(s,s)GS2xS2:d[s,S]>^} 
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is compact, and 

min \S-\s) - S-\S) \ =fii>0, 

(s,S)gP, 

Pi+1 C Pi,fii+i > fii,0 < i < h. Fix io = h - 2. Then |Ki(to) - i^iiti)\ > if 
i = io,io + 1. 

Similarly, we also have that there exists a positive number hq < 1 such that 
\i^i{to) - Ki(ti)\ < fiQ if < i < h. 

Without loss of generality assume that Ki(to) < i^i{ti) for all i. Thus, there is 
a number k among io, io + l such that Hk{ti) — Kk(to) > f^i^ and Kk(to) — Kk-i{ti) > 
i(l -/io). 

One can easily check that (C8) implies that z/i = — 2LMi > and 
z/2 < z/i, where z/2 = + ^V- 

We shall show that constants eo,ei for Definition 13.11 can be taken equal to 
eg = im?7, ei = min{/ii(,, ^(1 - /io)}- 

Assume that ||z(Cfc(to)) - ^i(Cfc(^o))|| < i^i- 

Then, fort G [^(to), Cfel^i)], 

= e^(*-C^(*"))(X+i?)z((a(to))+ r e^(*-^)/(;^(.))rfs+e-^(*-^'=(*«))H^(^((a(to))), 

^i(t) = e^(*-^'=(*»))zi((a(to)) + r e^(*-^)/(zi(s))f/s. 

-'Cfe(to) 

We have that 

Mt) - z,it)\\ = ||e^(*-^'=(*<'))[5^((a(to)) + WiziiCkitom + 

mrj — m{vi + 2LMi) > eo- 
If \\z{Ck{to)) - Zi{Ck{to))\\ > U2, then, for t G [Cfc_i(ti), Cfc(to)], 

-^Cfc(to) 
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and 

\\z{t) - zi{t)\\ > mz/2 - m2LMi = eo- 

The theorem is proved. 

On the basis of Theorems I3.m3.3l we can conclude that the following theorem 
is valid 

Theorem 3.4 // conditions (CI) - (C8) are fulfilled, then VCB is a chaotic 
attractor. 

Next, we are going to show that the attractiveness arranged by Theorem 12.21 
causes an amusing phenomena that the chaotic properties being attributed to all 
solutions of ([1]), not only those from VCB. 

3.2 The chaos of the initial value problem 

Let us first introduce the following definitions. 

Definition 3.5 We say that (QP is sensitive on A if there exist positive real 
numbers Kq,Ki such that for every solution z{t) = z(t,tQ, ZQ),to G A, o/ (QP and 
for each 6 > 0, H > 0, one could find a number ti G A, such that |to~^i| < and 
an interval Q from [tg, oo) with length no less than Ki such that \ \zi{t) — z{t)\ \ > 
K,o,t G Q, for all solutions Zi(t) = z(t,ti,zi), \ \zo — zi\\ < H, of and there 
are no points of discontinuity of zi{t) and z{t) on Q. 

Definition 3.6 The set of all periodic solutions (f){t,tQ),to g A, of (J\) is called 
dense in VC if for every solution z{t) = z{t,ti, zo),ti G A, of the initial value 
problem, and each e > 0, E > 0, there exists a periodic solution (pit) and an 
interval J C with length E such that 0(t)(e, J)z{t). 
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Definition 3.7 A solution z{t) = z{t,t* , zo),t* G A,t > t*, of HP is called 
dense in the set of all orbits of VC if for every solution zi{t) G PC, and each 
e > 0, -E > 0, there exists an interval J with length E and a real number ^ such 
that z{t + i){e, J)zi{t). 

Definition 3.8 Problem is chaotic on A if 

(i) it has the chaotic attractor VCB; 
(a) it is sensitive on VC; 

(Hi) the set of all periodic solutions (pit, to), to G A, is dense in VC; 

(iv) there exists a solution z{t,to, zo),to G A, of (Op, which is dense in the set 
of all orbits ofVC; 

Theorem 3.5 Assume that conditions (CI) — (C8) are fulfilled. Then (QP is 
chaotic on A. 

Proof: Condition [i] of Definition 13.81 is verified in Theorem 13. 4[ Using the 
attractiveness of VCB one can verify that parts {ii), {Hi) and {iv) are also vahd. 
Indeed, let us show first that {iv) is true. We will show that the solution z^ {t) = 
z{t,t*) dense in VCB by Theorem 13.21 is also dense in VC. Fix a solution z{t) = 
z{t,to,zo) G VC and a positive e. The solution z{t) is attracted by the bounded 
solution z{t,to) G VCB,z{t) ^ z{t,to), so that ([9]) is true. Hence, one can find a 
positive number T such that \\z{t,to, Zq) — z{t,to)\\ < e/2 if t > T. In the proof 
of Theorem 13.21 it was shown that there exists an arbitrarily large number Ei, 
where we may assume that Ei > 2T, and a natural number m such that z*(t + m) 
is e/2— equivalent to z{t, to) on the interval J = {Ei/2, Ei). Now, it is easy to see 
that z^..{t + m){e, J)z{t). So, {iv) is prowen. Part {Hi) can be verified in a similar 
way. Now consider condition {ii). To prove it, we will follow the proof and the 
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notation of Theorem 13.21 and we will show that the numbers in Definition 13.51 
can be as kq = eo/2 and Hi = ei, where cq and ei are the numbers used in the 
proof of Theorem 13. 2[ Consider a solution z(t,tQ,ZQ) G VC. From ([9]) it follows 
that for a fixed H > there exists a number T = T{H, zq) such that if t > T, 
then \\z{t,^,zi) — z{t,^)\\ < eo/4, where ^ G A is arbitrary, and \\zi — Zo\\ < H. 
Analyzing the proof of Theorem 13.21 one can readily observe that the number n, 
which satisfies all other demands of the proof, can be chosen greater than T + 2. 
Consequently, the number ti G A corresponding to h satisfies \ti — to I < ^) aiid 
\\z{t,tQ) -z(t,ti)|| > eo,t G Q. Hence, ||z(t, to, ^^o) - z{t,ti,Zi)\\ < ||z(t,to) - 
z{t,ti)\\ - \\z{t,tQ,ZQ) - zit,to)\\ - \\z{t,ti,zi) - z(t,ti)\\ > eo/2 = Ko,t e Q, 
where interval Q has length not less than Ki. The theorem is proved. 

Remark 3.1 In J^7| / the author considers the following initial value problem for 
a system of impulsive differential equations: 

x' = A{t)x + G,t=^ ri{m), 

x(0) = xo, (15) 

where t G [0, oo),a; G M", the vector xq is fixed, G and I are constant vectors, 
and 

(Dl) the associated homogeneous system x' = A{t)x is uniformly exponentially 
stable; 

(D2) an impulsive interval sequence Ti{m),i = 1,2,..., is defined by the formula 

i 

Ti = ^{rrik + 9), 

k=l 

where nik, k > 1, are iterations of a sensitive map h defined on the bounded 
interval H C M, such that there exists a number 6' > satisfying 6' + inf{t : 
t G S} > ^' > 0. 
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The following definition was given in the paper. 

Definition 3.9 141^ The solution of fT5\) is said to be chaotic with respect to the 
impulsive interval sequence provided that (i) this solution is uniformly bounded 
on the interval [0, cxd); (ii) there exist two positive numbers eo,cr > 0, such that 
for any m G S one can find a point m* G S arbitrarily close to m, and a positive 
number T satisfying 

/^({^o ^ t < T : \\x(t, m) — x{t, m*) \\ > a}) > cq, 
where /i(-) represents the Lebesgue measure of a given set. 

The author proves that under certain conditions, including (Dl), (-D2), the solu- 
tion of ( [73]) is chaotic in the sense of Definition \3. 9i 
Consider the following problem 

x'{t) = 0,t=^Ti{m), 
Ax|t=^^(„) = (-1)*, 

x{0) = xo, (16) 
where t,x eM, xq is fixed, and 

i 

Ti = ^^(mj + 1), mi = m,m E (0, 1), rrii = h{mi^i), h{ra) = 0. 
1=1 

It is easily seen that Ti = i + m — l,i = 1,2,... Consider two solutions of the 
problem, x{t,m),x{t,m*),m,m* G (0,1), m ^ m*. Obviously, yu({0 < t < k : 
\x{t, m) — x{t, m*)\ > 1} > k\m — m*\, where k is a positive integer, and solution 
x{t,m),m G (0,1), is bounded on [0,oo), as is x{t,m*). In other words, the 
solution is chaotic by Definition \3. 9\ even though there is no "chaotic" assumption 
on the map h, analogous to (-D2), and no condition similar to (Dl) is imposed on 
the system. So, this example shows that a more scrupulous investigation of the 
problem is required. In particular, not only sensitivity, but also other ingredients 
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of chaos should be involved in the discussion. We additionally emphasize that 
sensitivity in IJl^ is defined with respect to the impulsive interval sequence, or, it 
is better to say, with respect to the parameter m G S. Consequently, it cannot be 
considered a chaotic property. Our definition of sensitiveness involves the initial 
data, and it is natural for dynamics |^ [2B E^- Finally, we should remark that 
system ( f7g|) is not sensitive in the sense of Definition \3.5[ 

3.3 The period-doubling cascade 

Now we consider /i > 0, /i being the parameter for the logistic map. 

Since we use the logistic map as a tool to create the chaos, we expect to 
observe the period-doubling cascade. 

Consider the logistic map h{t., fi) = — t), assuming this time that /i < 4. 
It is known [10], that there exists an infinite sequence 3 < yUi < < • • • < 
/ifc . . . < 3.8284. . . such that h{t,fii),i > 1, has an asymptotically stable prime 
period— 2* point t* with a region of attraction (t* — 6i,t* + 6i). And beyond the 
value 3.8284. . . , there are cycles with every integer period [9]. 

By the results of our paper the period-doubling process is occurring in the 
bounded region ||a;|| < Mi of the space M", common for all /i > 0. The chaotic 
properties are observed in the region, and, by Theorem 12.21 bounded solutions 
from this domain attract all other solutions. Moreover, the cascade generates 
infinitely many periodic solutions. 

Example 3.1 Consider the following initial value problem 

x[ = 2/5x2 + ^ sin^ X2, 

x'2 = 2/5xi + /sin^a;i,t 7^ Ci(^o), 

Axi\t=(^^to) = -32^1. 

Ax2\t=Q(to) = -lx2 + W{x2), (17) 
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where W{s) = 1 + s^, if \s\ < 1,1 is a positive constant, and W{s) = 1 + if 
\s\ > I. One can easily see that all the functions are lipschitzian with a constant 
proportional to I. The matrices of coefficients are 

A=i «J/5y -4/3 



2/5 y ' V - 4/3 

The matrices commute, and the eigenvalues of the matrix 



A + Ln{l + B) 



-ln3 2/5 
2/5 - In 3 



are negative: Ai^2 = — ln3 ± 2/5 < 0. 

Let us check if condition (C3) holds. It is sufficient to verify that — + 
W{s) ^ 0, for all s G M. 

If \s\ < I, then — |s + W{s) = — |s + 1, and it is never equal to zero. 

If \s\ < I, then — |s + W{s) = P — ^x + 1. For the last expression to be zero, 
we need, s = —2(1 + /^), that is \s\ > I. We have a contradiction. Thus, condition 
(C3) is valid. All the other conditions required by the Theorems could be easily 
checked with sufficiently small coefficient I. That is for /x > 4 there is the chaotic 
behavior. 

The numerical simulation of the chaos is not an easy task since even the 
verification of sensitivity requires two close values of the initial moment in the 
Cantor set A, which cannot be found easily. Hammel et all. fj^ have given a 
computer- assisted proof that an approximate trajectory of the logistic map can be 
shadowed by a true trajectory for a long time. This result and the continuous 
dependence of the solutions on the sequence of discontinuity points make possible 
the following appropriate simulations. 

In what follows we will demonstrate for /x = 4 and fi = 3.8282 the sensitivity 
and intermittency respectively. 

Assume that I = 10~^. Choose fi = 4 in [Tl\) and consider two solutions 



x{t) = {xi,X2),x{t) = (xi,X2), with initial moments to = 7/9 and to = 7/9+3 
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respectively. That is, we take the initial values close to each other, and, moreover, 
the solutions with identical initial values, x{to) = x(to) = (0.005,0.002). The 
graphs of the coordinates of these solutions (Fig. 1) show that the solutions 
abruptly become different when t is between 15 and 20, despite being very close to 
each other for allt in the interval (to, 15). One can conclude that the phenomenon 
of sensitivity is numerically observable. 

Moreover, if one consider the sequence {xi{n),X2{n)),n = 1, 2, 3, ... , 75000, 
in Xi,X2— plane (Fig. 2), then the chaotic attractor can be seen. 



/ / , / / , I II i I !l I li I I I i ! II II I j 1 1 

Uliihllil II I I lillllljl ^1 

III'! r I n'i'i , II / ir'i'^r II' ii I'' 



^r'lUI'il iil^l'i 



time t 
(a) 



0.5 1 — 
0.45 - 

0.4 - 
' 0.35 - 

0.3 - 
0.25 - 

0.2 — 



time t 
(t>) 



Figure 1: Simulation results. Graphs of coordinates Xi(t),X2(t) are in the blue, 
and of coordinates Xi{t).,X2{t) are red. The coordinates in (a) and (6) abruptly 
become significantly different when t is near 15, while almost coinciding for all t 
in the interval (ti, 15). 



Consider n = 3.8282 in p7| j. The coefficient's value o/ 3.8282 is such that 
the logistic map admits intermittency fj^ . Fix to = 0.5 and take a solution 
{xi{t),X2{t)) of the last system with the initial condition xi(to) = 0.002, ^2(^0) = 
0.005. The result of simulation can be seen in Fig. 2. 
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Figure 2: The chaotic attractor by a stroboscopic sequence X2(^)), 1 < 

n < 75000, is observable. 
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Figure 3: Simulation results. In (a) and (6) for coordinates Xi and X2 respectively 
the intermittency phenomenon is indicated. 
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4 Conclusion 



The complex dynamics is obtained using Devaney's definition for guidance. We 
also prove the existence of a chaotic attractor consisting of solutions bounded on 
the whole real axis, and observe the period-doubling cascade for the problem. 
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